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The nuclear-mass dependence of azimuthal cross section asymmetries with respect to charge and 
longitudinal polarization of the lepton beam is studied for hard exclusive electroproduction of real 
photons. The observed beam-charge and beam-helicity asymmetries are attributed to the interfer- 
ence between the Bethe-Heitler and deeply virtual Compton scattering processes. For various nuclei, 
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the asymmetries are extracted for both coherent and incoherent-enriched regions, which involve dif- 
ferent (combinations of) generalized parton distributions. For both regions, the asymmetries are 
compared to those for a free proton, and no nuclear-mass dependence is found. 

PACS numbers: 13.60.-r, 13.60.Le, 13.85.Lg, 14.20.Dh, 14.40.Aq 



Lepton scattering experiments constitute an impor- 
tant source of information for understanding nucleon 
structure in the context of QCD. Until recently, this 
structure was described by two categories of non- 
perturbative objects, form factors and parton distribu- 
tion functions (PDFs), which have been measured in elas- 
tic and deep-inelastic scattering (DIS) experiments, re- 
spectively. In the last decade, generalized parton dis- 
tributions (GPDs) have been recognized as a key 
concept for the description of hard exclusive processes. 
GPDs offer a multi-dimensional representation of the 
structure of hadrons at the partonic level, correlating the 
longitudinal momentum fraction carried by the parton 
with its transverse spatial coordinate 0-0] ■ For recent 
theoretical reviews, see Refs. [IMl- 

Generalized parton distributions depend on the 
squared four-momentum transfer t to the nucleon and on 
x and £, which represent respectively the average and half 
the difference of the longitudinal momentum fractions 
carried by the probed parton in initial and final states. 
Nucleon elastic form factors and parton distribution func- 
tions appear as x-moments and kinematic limits (for 
t, £ — > 0) of GPDs, respectively. The skcwness param- 
eter £ is related to the Bjorken variable xb = Q 2 / (2M v), 
as £ f=a xb/(2 — x-q) in the Bjorken limit where Q 2 — > oo at 
fixed values of Xb and t. Here, M is the target mass and 
— Q 2 is the squared four-momentum of the exchanged 
virtual photon with energy v in the target rest frame. 
Most often discussed are the four twist-2 quark-helicity- 
conserving GPDs for each quark species in the nucleon: 
the quark-polarization averaged distributions H and E 
and the quark-polarization related distributions H and 
E. 

Among all presently practical hard exclusive probes, 
deeply virtual Compton scattering (DVCS), i.e., the hard 
exclusive leptoproduction of a real photon, appears to 
have the most reliable interpretation in terms of GPDs. 
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The final state of the DVCS process in which the real 
photon is radiated by a quark is intrinsically indistin- 
guishable from that of the Bethe-Heitler (BH) process in 
which a real photon is radiated by the incoming or outgo- 
ing lepton. Access to the DVCS amplitude is provided by 
interference between the Bethe-Heitler and DVCS pro- 
cesses, e.g., via the measurement of the cross-section 
asymmetries with respect to the lepton beam helicity and 
charge. 

This paper reports the first experimental study of 
DVCS on nuclear targets. Nuclei provide a laboratory 
where, compared to the free nucleon, additional informa- 
tion can be obtained on GPDs by observing how they 
become modified in the nuclear environment. Therefore, 
studies of nuclear GPDs offer a new opportunity to in- 
vestigate the nature of the nuclear environment. 

In lepton-nucleus scattering, two processes can be dis- 
tinguished for both DVCS and BH: (a) the coherent pro- 
cess where the electron scatters off the whole nucleus, 
which stays intact; (b) the incoherent process where the 
electron scatters quasi-elastically from an individual nu- 
cleon, breaking up the nucleus. 

For coherent scattering, various DVCS observables 
have been estimated theoretically [3, EH- In these es- 
timates, nuclear GPDs are expressed in terms of nucleon 
GPDs convoluted with the distribution of nucleons in the 
nucleus. The t dependence is modeled using nuclear elas- 
tic form factors. These models predict an enhancement 
of the beam-charge and beam-helicity asymmetries for 
spin-0 and spin-1/2 nuclei compared to the case of a 
free proton. 

Recently, coherent DVCS on nuclei has been suggested 
to provide new insights into the origin of the EMC ef- 
fect as models that attempt to explain the EMC 
effect in the forward case (t, £ — > 0) also predict nuclear 
GPDs that differ from those of a free nucleon ('general- 
ized' EMC effect). GPD models embodying PDFs that 
describe the EMC effect observed in inclusive DIS pre- 
dict a much larger generalized EMC effect for DVCS ob- 



servables Ulj]. In Rcf. 22 1, this enhancement is at 



tributed to the transverse motion of quarks in nuclear tar- 
gets, while Ref. 23[ relates the enhancement to mesonic 
degrees of freedom in hard reactions on nuclei, which have 
been invoked in the 'pion excess' models to explain the 
EMC effect in inclusive DIS 0, El 0] . An observable 
found to be sensitive to mesonic degrees of freedom is the 
real part of the DVCS amplitude, which is predicted to 
strongly depend on the nuclear mass number A 23 1. 

Incoherent scattering is approximated by scattering 
on free nucleons. In the kinematic conditions of this 
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experiment, scattering on the proton dominates due to 
the fact that the BH process dominates the single photon 
production rate and the BH process on the neutron is 
suppressed because of the small electromagnetic form 
factors compared to those of the proton. Therefore 
the asymmetries for nuclei in the incoherent case are 
anticipated to be similar to those for the proton. The 
role of the neutron contribution was studied in Ref. [25| . 
It was shown to decrease the asymmetries measured in 
incoherent nuclear DVCS at larger values of —t. 

The cross section for hard exclusive leptoproduc- 
tion of real photons reads 



da 



x B e 



ITI 



dx B dQ 2 d\t\ 



32(2^) 4 Q 4 vTT?' 



(1) 



where e represents the elementary charge, e = 2x B M/Q 
and T is the total reaction amplitude. The azimuthal 
angle </> is defined as the angle between the lepton scat- 
tering plane and the photon production plane spanned by 
the trajectories of the virtual and real photons, following 
Ref. [26j. The scattering amplitudes of the DVCS and BH 
processes add coherently. The cross section is then pro- 
portional to the squared photon-production amplitude 
written as 



|T| 2 = |T DV cs| 2 + |Tbh| 2 +L 



where the interference term I is given by 



I = T 



dvcsTbh 
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DVCSTBH 
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The BH amplitude Tbh is calculable from measured clas- 
tic form factors of the (nucleon) nucleus when modelling 
the observablcs for the (in)coherent process. At leading 
order in the fine structure constant a and for an unpo- 
larized target, the squared BH amplitude |Tbh| 2 is in- 
dependent of beam polarization and lepton charge. In 
contrast, the squared DVCS amplitude |Tdvcs| 2 & n d the 
interference term I depend on the beam helicity, while the 
interference term depends also on the lepton charge. For 
a longitudinally polarized lepton beam and unpolarized 
target, these dependences read fu\ 



BH 



- DVCS I 



cosine 
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(5) 



J2 K cos ( r 



n=0 



P b [s\ sin 4> + 4 sin(20)] ) . (6) 



Here, Pb denotes the longitudinal beam polarization, eg 
the beam charge in units of the elementary charge, V\ (4>) 
and V2{4>) are the known ^-dependent lepton propaga- 
tors in the BH process, and the kinematic factors read 
Kbh = l/[x B t(l + e 2 ) 2 ] and Kj = l/{x B yt) with y the 
fraction of the incident lepton energy carried by the vir- 
tual photon in the target rest frame. The dependences of 
the coefficients c„ and s n on GPDs are given in Ref. ll| 1 
for a spin-1/2 target and in Ref. [27| for a spin-0 target. 
For a spin-1/2 target, and within the typical kinematic 
conditions of this experiment, the coefficients related to 
only twist-2 quark GPDs appearing in the interference 
term can be approximated as 



c\ cx PiRe-H, 



s\ cx FiIm'H, 



(7) 

(8) 
(9) 



where H denotes the Compton form factor that is a 
convolution of the GPD H with the hard scattering 
amplitude, and Pi is the Dirac form factor. 

In this paper we present a study of hard ex- 
clusive production of real photons in the reaction 
eA -> e^X. The data were collected with the HERMES 
spectrometer during the period 1997-2005. The 
27.6 GeV HERA electron or positron beam at DESY 
was scattered off gaseous hydrogen, helium, nitrogen, 
neon, krypton, and xenon targets (see Table [TJ) . (Results 
from a deuterium target will be reported elsewhere (29|.) 
The HERA beam was transversely self-polarized due 
to the Sokolov-Ternov mechanism [30] • Longitudinal 
polarization of the beam was obtained by using a pair 
of spin rotators located before and after the interaction 
region of HERMES. The beam helicity was reversed 
every few months. The beam polarization was measured 
by two independent HERA polarimeters [HI, [32| with 
a combined fractional systematic uncertainty of up to 
3.4%. This analysis makes use of the full data set with 
nuclear targets and a subset of data with a hydrogen 
target taken in the years 2000 and 2005 corresponding 
to approximately 130 pb -1 (100 pb -1 ) for the positron 
(electron) sample. (The results from the full 1996-2005 
hydrogen data set has been reported elsewhere [HI].) 
For hydrogen, krypton and xenon targets, data for both 
positron and electron beams are available. 

A brief description of the event selection is given here. 
More details can be found in Refs. 13, 3f|. Events were 
selected if exactly one photon and one charged track iden- 
tified as the scattered lepton were detected. The hadron 



1 Note that the azimuthal angle 
the one used in Ref. Illl (<f> = it 



defined here is different from 
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1/2 


227 


0.50 


-0.51 


He 


n 

yj 




U.JU 




N 


1 


51 


0.39 


-0.40 


Ne 





86 


0.52 


-0.55 


Kr 





77 


0.43 


-0.41 


Xe 


0, 1/2, 3/2 


47 


0.32 


-0.38 



TABLE I: Targets used for this analysis, their spins, the 
corresponding integrated luminosity L, and the average po- 
larization for the two helicity states of the beam. Note that 
the xenon target is composed mainly of the isotopes 129 Xe 
(spin-1/2), 131 Xe (spin-3/2) and 132 ' 134 Xe (spin-0) with frac- 
tional contributions of 26%, 21% and 36%, respectively. For 
all other targets, the admixture of isotopes with spin different 
from that given in the table is less than 10%. 



Monte Carlo sum 
coherent BH 
incoherent BH 
associated BH 
semi-inclusive 




10 12 14 

M£ [GeV 2 ] 



contamination in the lepton sample is kept below 1% by 
combining the information from a transition-radiation 
detector, a prcshowcr scintillator detector, and an elec- 
tromagnetic calorimeter. The kinematic requirements 
imposed are 1 GeV 2 < Q 2 < 10 GeV 2 , 0.03 < x B < 
0.35, v < 22 GeV, and W > 3 GeV, where W is the in- 
variant mass of the virtual-photon nucleon system. The 
real photon is identified by a 'neutral cluster', which is 
defined as an energy deposition larger than 5 GeV in the 
calorimeter and larger than 1 MeV in the preshower de- 
tector, and the absence of a corresponding charged track. 
The angular separation (9 7 » 7 between the virtual and real 
photons is required to be larger than 2 mrad. This value 
is chosen in order to optimize the combined systematic 
and statistical uncertainties for the asymmetries due to 
the degraded <p resolution at low 7 * 7 and the enhanced 
production of real photons on nuclear targets in the small 
(9 7 * 7 region 2 [3a |. An upper bound of 45 mrad is im- 
posed on this angle in order to improve the signal-to- 
background ratio. 

The recoiling system was not detected. Instead, an 
'exclusive' sample of events is selected by requiring the 
squared missing mass = (q + p — q') 2 to correspond 
within the experimental resolution to the squared proton 
mass. Here, q (q 1 ) is the four-momentum of the virtual 
(real) photon and p = (M p , ) with M p the proton mass. 
This selection criterion is chosen by means of a Monte 
Carlo (MC) simulation of the missing mass distribution. 
The result of the simulation is shown in comparison with 
the experimental data in Fig. [T] In the MC simulation 



the expressions in Eqs. 35 and 36 of Rcf. [llj are used 
for the incoherent BH process. The simulation also takes 
into account the incoherent BH process where a nucleon 
is excited to a resonant state (known as associated pro- 



FIG. 1: (color online) Distributions in squared missing mass 
from data using positron (filled circles) or electron (empty 
circles) beams and a xenon target compared to a MC simula- 
tion (solid line). The latter includes coherent Bethe-Heitler 
(BH) (dashed line), incoherent BH (short-dashed line) and 
associated BH (filled area) processes as well as semi-inclusive 
background (dash-dotted line). The two vertical solid lines 
enclose the selected exclusive region for the positron data. See 
text for details. 



duction) using a parameterization of the total "f*p cross 
section for the resonance region from Rcf. 



381 and cal- 



culating the individual cross sections for single-meson 
decay channels, e.g., A+ -> pn , with the MAID2000 
program [39j]. For the coherent BH process, the param- 
eterizations of the form factor for the respective nuclear 
targets are taken from Refs. [2|| 37 1. In addition, semi- 



Note that this value is the only difference from earlier HERMES 



analyses, for which 



> 5 mrad. 



inclusive production of neutral mesons (mostly 7r°) is in- 
cluded, where either only one photon from the 7r° — > 77 
decay is detected or these photons cannot be experi- 
mentally resolved. For this process, the MC generator 
LEPTO [Io[ is used in conjunction with a set of JET- 
SET [4l[ fragmentation parameters that had previously 
been adjusted to reproduce multiplicity distributions ob- 
served by HERMES 0. Not included in the simula- 
tion is radiation of more than one photon, which would 
move events from the peak to the continuum, nor the 
DVCS process. The latter contribution is highly model- 
dependent. In the GPD model used in Ref. [43j it varies 
between 10% and 25% of the BH yield for production 
from a hydrogen target (3) ■ 

The 'exclusive region' for the positron data is defined 
as -(1.5 GeV) 2 < M 2 < (1.7 GeV) 2 , where the lower 
limit is chosen to be three times the resolution in M\ 
from the squared proton mass, and the upper limit to 
optimize the signal-to-background ratio. Since the M\ 
spectrum of the electron data is found to be shifted by 
approximately 0.18 GeV 2 towards smaller values relative 
to that of the positron data, the exclusive region for elec- 
tron data is shifted accordingly. One quarter of the effect 
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coherent BH 

ncoherent BH 
associated BH 




0.1 0.12 

-t [ GeV 2 ] 



FIG. 2: (color online) Distribution (points) of events selected 
in the exclusive region as function of —t compared to a MC 
simulation (solid line). The latter includes coherent Bethe- 
Heitler (BH) (dashed line), incoherent BH (dotted line) and 
associated BH (filled area) processes. Background from semi- 
inclusive neutral meson production is not included. 



of this shift on the results presented below is assigned as 
a contribution to the systematic uncertainty. 

As the recoiling system was not detected, t is inferred 
from the measurement of the other final-state particles. 
For elastic events, the kinematic relationship between the 
energy and direction of the real photon permits the cal- 
culation of t without using the measured energy of the 
real photon, which is the quantity subject to larger un- 
certainty. Thus the value of t is calculated as 



t = -Q 2 -1v(v- \Jv 2 + Q 2 cos e Tl ) ^ 

1 + ik; ( v ~ v 7 ^ 2 + Q 2 cos6> 7 . 7 ) 

for the exclusive event sample. The error caused by ap- 
plying this expression to inelastic events is accounted 
for in the MC simulation that is used to calculate the 
fractional contribution of background processes per kine- 
matic bin. The quantity —t is required to be smaller than 
0.7 GeV 2 . 

Coherent scattering on nuclear targets is separated 
from incoherent scattering by exploiting its characteristic 
t dependence. For both DVCS and BH, coherent scatter- 
ing occurs at small values of —t and rapidly diminishes 
with increasing \t\. However, a complete separation of 
the two scattering processes is impossible at HERMES. 
Cohcrent-enric/ierf and incoherent-enric/ied samples are 
selected according to a —t threshold that is chosen to vary 
with the target such that for each sample approximately 
the same average kinematic conditions are obtained for 
all targets. The kinematic distributions of elastic coher- 
ent and incoherent processes are determined using the 
MC simulation described above and presented in Figs. Q] 
and [2] for xenon, as an example. The t distribution of 
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*± 
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-t < 0.043 


-0.018 


0.068 


1.73 


66 
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Ne 


-t < 0.044 


-0.018 


0.068 


1.74 


68 
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Kr 


-t < 0.070 


-0.018 


0.064 


1.63 


69 
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Xe 


-t < 0.078 


-0.018 


0.062 


1.60 


66 


4 



TABLE II: Average kinematics and fractional contributions 
from coherent processes (purity) and associated processes in 
the coherent-enriched sample for the various targets. 
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t threshold 
[GeV 2 ] 


(t) 
[GeV 2 ] 


(xb) 


(Q' J ) 

[GeV 2 ] 


% of 
incoh 
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assoc. 


H 


-t > 0.077 


-0.200 


0.109 


2.89 




20 


He 


-t > 0.084 


-0.200 


0.107 


2.78 


61 


28 


N 


-t > 0.083 


-0.200 


0.113 


2.93 


60 


28 


Nc 


-t > 0.075 


-0.200 


0.111 


2.92 


65 


28 


Kr 


-t > 0.067 


-0.200 


0.108 


2.84 


57 


30 


Xe 


-t > 0.060 


-0.200 


0.107 


2.86 


56 


30 



TABLE III: Average kinematics and fractional contributions 
from incoherent processes (purity) and associated processes 
in the incoherent-enriched sample for the various targets. 



events selected in the exclusive region is shown in Fig. [5] 
for xenon together with the simulated contributions of 
coherent and incoherent processes. Tables |H] and IIIII 
summarize the average kinematic conditions for the var- 
ious targets for the coherent-enriched and incoherent- 
enriched samples, respectively, and give their purities de- 
fined as fractions of the total simulated yield. Also shown 
for each sample is the simulated fractional contribution 
from the associated BH process. For hydrogen, kryp- 
ton and xenon, the coherent-enriched region is further 
explored as a function of t (see Table ITV)) . 

As for these data coherent scattering could not be iden- 
tified event-by-event, kinematic variables that depend 
on the target mass are always calculated using the proton 
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t range 


(t) 
[GeV 2 ] 


(xb) 


(Q 2 ) 

[GeV 2 ] 


% of 
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[GeV 2 ] 




coh 


assoc. 
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0.0 <-t< 0.008 


-0.006 


0.054 


1.38 




< 1 


Kr 


0.0 <-t< 0.010 


-0.006 


0.053 


1.37 


92 


< 1 


Xe 


0.0 <-t< 0.010 


-0.006 


0.053 


1.37 


92 


< 1 


H 


0.008 <-t< 0.020 


-0.014 


0.069 


1.75 




1 


Kr 


0.010 < -t < 0.020 


-0.014 


0.064 


1.63 


71 


2 


Xe 


0.010 <-t< 0.020 


-0.014 


0.062 


1.67 


71 


2 



TABLE IV: Average kinematics and fractional contributions 
from coherent processes (purity) and associated processes in 
two t subranges of the coherent-enriched sample for hydrogen, 
krypton and xenon. 
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mass. This does not influence the selection of exclusive 
events since the values of the relevant kinematic variables 
calculated using the proton mass are highly correlated 
with those calculated using the actual target mass. Also, 
the calculation of the t value is affected negligibly. 

The full cross section for exclusive production of real 
photons on unpolarizcd targets (U) by a longitudinally 
polarized beam (L) can be written as 

a(4>) = cruu,o(</>)[ 1 + eeAc (</>)+ -Pb-4,LU,Dvcs(0) 

+ e,P b -4 L u,i (</>)], (11) 



where ctuu,o(</>) = 



is the 



cross section for an unpolarizcd target averaged over both 
beam charges (+, — ) and over both positive (— >) and neg- 
ative (<— ) beam helicities. The beam-charge asymmetry 
Ac and beam-helicity asymmetries *Alu,dvcs and *Alu,i 
are defined in Eqs. Q~2l [14J and \W{ respectively. Each 
definition is complemented by the corresponding relation 
to the coefficients given in Eqs.Q][G] 



A C (<P) = 



a ' 



(12) 



<nju,oW>) V l {4>)V 2 W>)^ 



- (T ^~ 
3 

^ c l n cos(n0),(13) 



A 



LU,DVCS 



(0) 



1 



DVCS 



O"UU,o(0) 



sm ( 



— a 



+ (J- 



(14) 
(15) 

(16) 



h o— >- 

2 

^ s\ sin(n0X.17) 



As the term of Eq. QTJ including -4lu,i depends on both 
beam hclicity and beam charge, the DVCS and interfer- 
ence beam-helicity as ym metries can be separated. Such 
a combined analysis [451 ] was performed for hydrogen, 
krypton and xenon, where data for both electron and 
positron beams are available. The asymmetries defined 
in Eqs. [321 El and [16] are expanded in terms of the 
following harmonics in <fi: 



AcW 
-4lu,dvcs(0) 



-4 



s (04>) 



(18) 



A c ° s{2 * ] cos(2$ + A c ° s{3<p) cos(3^), 



^LU.DVCS bm ' 



Ai«fhin(2<f>)- 



(19) 
(20) 



Using the method of maximum likelihood, the Fourier 
coefficients A, hereafter called asymmetry amplitudes, 
are simultaneously extracted from the event yield that 



is proportional to the cross section of Eq. QT] Although 
these asymmetry amplitudes differ somewhat from the 
coefficients in Eqs. Q31 HSl andfTTl they are well defined 
and can be computed in various GPD models for direct 
comparison with data. 

For helium, nitrogen and neon, only data with 
a positron beam were collected. The single-charge 
(positron) beam-helicity asymmetry is defined as 



•Alu.+ M 



(7 + (J 
1 



1 



C r UU,+ (0) 

-1 



s ? vcs sin< 



e e Ki 



(21) 
(22) 



XX sin(n0), 



where <7uu,+ (0) = \{ a ^ + < ji ~)- I 11 this case the event 
yield that is proportional to the cross section of Eq. [11] 
is fitted by 



A 



LU, 



A LU + Sin< 



.4 : 



sin(2^> 
LU,+ 



sin(20). (23) 



This method does not allow for a separation of squared 
DVCS amplitude and interference term in the beam- 
helicity asymmetry It was used in an earlier extraction 
of beam-helicity asymmetries for hydrogen [46^ . 

In each kinematic bin, the extracted asymmetry am- 
plitudes are corrected for background from the decay of 
semi-inclusivcly produced neutral mesons, mainly pions. 
The corrected asymmetry amplitude is then obtained as 



A c 



^■raw /semi ' -^s 
1 /semi 



(24) 



where A raw stands for the extracted raw asymmetry am- 
plitude and /semi and A sem \ for the fractional contribution 
and corresponding asymmetry amplitude of the semi- 
inclusive background, respectively The background con- 
tribution /semi, estimated from MC simulations, ranges 
from 1% to 11% depending on the kinematic conditions 
and amounts to 3.5% on average. Since the beam- 
chargc-dcpcndcnt background asymmetry is zero at lead- 
ing order QED, the semi-inclusive background consti- 
tutes a dilution for Ac and effectively also for .Alu.Ij 
since it cancels in the latter case. In order to correct 
-4lu,dvcs and Altj for the semi-inclusive background, 
the size of the corresponding beam-helicity asymmetry 
is extracted from data by reconstructing neutral pions 
with a large fractional energy z = E^a/v > 0.8, as ac- 
cording to MC simulations only these contribute to the 
exclusive region [44j . These simulations show that the ex- 
tracted 7T° asymmetry does not depend on whether only 
one or both photons are in the acceptance. One half 
of the size of the full background correction is assigned 
as systematic uncertainty. Contributions from exclusive 
7T° production were found to be negligible at HERMES 
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FIG. 3: The cos <j> amplitude of the beam-charge asymmetry 
for hydrogen, krypton and xenon as function of t. The error 
bars (bands) represent the statistical (systematic) uncertain- 
ties. 



in a MC simulation based on a GPD model for exclu- 
sive meson production 1471 as well as in a data search for 
exclusive 7r° production [48| . Hence this conceivable con- 
tribution is not included in the systematic uncertainty. 

The asymmetry amplitude ^4 C orr, determined by ap- 
plying Eq. [24l is expected to originate from only elastic 
and associated production. Because essentially nothing 
is known about the asymmetry for associated production, 
no correction is made or uncertainty is assigned for the 
latter, but instead associated production is considered to 
be part of the signal in this analysis. The fractional con- 
tribution of associated processes is strongly t dependent, 
ranging from 3% in the lowest t bin to 50% in the highest 
t bin, with little dependence on A. 

Effects from detector acceptance, kinematic smearing, 
finite bin width, and from possible detector misalignment 
are estimated using a MC simulation based on the GPD 
model of Ref . [!|| . Note that a mistake has been found 
in this GPD model jE(J; however the model describes 
well the magnitude and kinematic dependences of previ- 
ously reported HERMES beam-charge [Hf and prelimi- 
nary beam-helicity asymmetries (5ll | and thus is consid- 
ered to be suitable for systematic studies. For each bin 
in —t, the asymmetry amplitude for hydrogen is (i) cal- 
culated at the mean kinematic values of a given bin and 
(ii) extracted from the reconstructed MC simulation that 
includes all experimental effects. The difference between 
these two amplitude values is included in the system- 
atic uncertainty. This uncertainty estimated for hydro- 
gen is applied to all targets. The validity of this approach 
was checked using MC simulations based on the model 
of Ref. that also parameterizes nuclear GPDs. The 
systematic uncertainty obtained for the nuclear targets 
is, within its statistical uncertainty, in good agreement 
with that estimated for hydrogen. 

The total systematic uncertainty is dominated by the 
uncertainty from the combined contributions of detec- 
tor acceptance, kinematic smearing, finite bin width, and 



FIG. 4: The sin0 amplitude of the beam-helicity asymme- 
try sensitive to the squared DVCS amplitude for hydrogen, 
krypton and xenon as function of t. The error bars (bands) 
represent the statistical (systematic) uncertainties. This am- 
plitude is subject to an additional 3.4% maximal scale uncer- 
tainty arising from beam polarimetry. 



from possible detector misalignment. This combination 
is added in quadrature with contributions arising from 
the background correction and the relative energy shift 
between the M\ spectra of positron and electron data. 
A scale uncertainty of up to 3.4% arising from beam po- 
larimetry is not included in the systematic uncertainty 
for the beam helicity related asymmetries. Also not in- 
cluded is any contribution due to additional QED ver- 
tices, as the most significant of these was estimated to 
be negligible in the case of helicity asymmetries [52j . 

In Figs. [3H5J the amplitudes of the beam-charge and 
beam-helicity asymmetries, A™ s< ^, -^ludvcS' -^lu I' an< ^ 
-^lu +> are s h° wn as functions of — t for unseparated co- 
herent and incoherent production. For the nuclear tar- 
gets, all other amplitudes in Eqs. [TBH2CII 1231 are found to 
be consistent with zero within 1.5 sigma of the statistical 
uncertainty. These other asymmetry amplitudes relate to 
coefficients that either embody higher twist quark GPDs 
or are kinematically suppressed as, e.g., the amplitude 
presented in Eq. [5] 

Figure [3] shows the amplitude A^? s for hydrogen, 
krypton, and xenon. The values for hydrogen from 
this a naly sis are consistent with those extracted previ- 
ously [111 |45[. For hydrogen, krypton and xenon, the 
availability of data with both beam charges allows for 
the separation of the azimuthal harmonics of the squared 
DVCS amplitude and the interference term. The beam- 
helicity amplitude -^ludvcS' snown m Fig. 01 is consis- 
tent with zero for all three targets over the full —t range. 
This is in agreement with the expected suppression of the 
amplitude. 

The beam-helicity amplitudes ^luI ano - -^lu+ snown 
in Fig.0 are substantial for all targets. For helium, nitro- 
gen and neon, where only positron beam data are avail- 
able, this amplitude receives also contributions from the 
squared-DVCS term. However, as the latter amplitude 
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FIG. 5: The t dependence of the sin <f) amplitude of the beam- 
helicity asymmetry sensitive to the interference term, ^lu'i, 
for hydrogen, krypton and xenon (full symbols) or to a linear 
combination of the interference and the squared DVCS ampli- 
tude, for helium, nitrogen and neon (open symbols). 
The error bars (bands) represent the statistical (systematic) 
uncertainties. This amplitude is subject to an additional 3.4% 
maximal scale uncertainty arising from beam polarimetry. 



is expected to be suppressed and found to be so for other 
targets in Fig. @J its contribution is assumed to be small 
here. 

The nuclear-mass dependence of the azimuthal beam- 
charge and beam-helicity asymmetries is presented sepa- 
rately for the coherent and incoherent-enriched samples 
in Figs. |B]and [7] The cos 4> amplitude of the beam-charge 
asymmetry is consistent with zero for the coherent- 
enriched samples for all three targets, while it is about 
0.1 for the incoherent-enriched samples without show- 
ing any dependence on the target mass within uncertain- 
ties. The sin <j) amplitude of the beam-helicity asymme- 
try shown in Fig.[7]has values of about —0.2 for both the 
coherent and incoherent-enriched samples without show- 
ing any dependence on A within uncertainties. In order 
to quantify nuclear effects, the asymmetry amplitudes 
for nuclear targets are compared to those for a free pro- 
ton. The ratio R\jj = ^lu^i +) a/^lu'i h °f the nuclear- 
to-hydrogen beam-helicity asymmetry amplitudes av- 
eraged over all targets, is found to be 0.91 ± 0.19 for 
the coherent-enriched sample and 0.93 ± 0.23 for the 
incoherent-enriched sample, both of which are compati- 
ble with unity. 

For incoherent scattering, the asymmetry for nuclei is 



expected to be similar to that for hydrogen aside from 
effects of the nuclear environment, since scattering on a 
proton dominates. Neglecting the neutron contribution, 
the value of i?LU for incoherent scattering is expected to 
be unity jlj|. In Ref. 25 1, the neutron contribution to 



incoherent nuclear DVCS is taken into account and i?LU 
is predicted to be between 0.66 and 0.74 at t = —0.2 
GeV 2 . Within the experimental uncertainties, the mea- 
sured ratio i?LU = 0.93 ± 0.23 agrees with both the ex- 
pected suppression of the neutron contribution in inco- 
herent scattering on nuclei and with the prediction of 
Ref. " 



25 



The results for the coherent-enriched samples can be 
compared to predictions based on simple models for nu- 
clear GPDs that express them in terms of nucleon GPDs 
14, [l5|. Within this approach, nuclear beam-charge 



and beam-helicity asymmetries are predicted to be es- 
sentially independent of A for heavier nuclei. Compared 
to the free proton asymmetry, nuclear beam-charge and 
beam-helicity asymmetries are expected to be enhanced 
for spin-0 and spin-1/2 nuclei. This predicted enhance- 
ment is based on the model-independent observation 
that DVCS takes place either on a proton or neutron in 
the nuclear target, while BH occurs predominantly only 
on a proton. The ratio i?LU of the nuclear-to-hydrogcn 
beam-helicity asymmetry amplitudes has been estimated 
in Ref. [14| for the pure coherent process to be about 



5/3 for spin-0 and spin-1/2 nuclei with Z=N, essentially 
independent of A. This value arises from the ratio of 
squared charges for an isoscalar to an isodublet state 
and the observation that for the valence quark PDFs 
d/u =1/2 in the kinematics of this experiment. For 
spin-1 nuclei, i?LU is predicted to be unity. Ref. Q 
also formulates a GPD model. Considering only leading 
twist GPDs and valence quark contributions, the pre- 
dicted value of A^u"^ for hydrogen is —0.26 for the kine- 
matic condition t = -0.2 GeV 2 , Q 2 = 2.5 GeV 2 and 
ib = 0.12. Including sea quark contributions, twist-3 
corrections and varying the main model parameters, the 
predicted amplitude is in the range 0.16 < |Aluj| < 0.37. 
In Ref. [l5| , a somewhat more elaborated calculation is 
presented where nuclear GPDs are expressed in terms of 
nucleon GPDs convoluted with the distribution of nucle- 
ons in the nucleus thereby accounting for nuclear binding. 
Within this approach, the ratio i?LU is predicted to be 
about 1.8 for neon and krypton for the kinematic condi- 
tion t = -0.018 GeV 2 , Q 2 = 1.58 GeV 2 and x B = 0.10. 

The nuclear beam-helicity amplitudes shown in Fig. [7] 
(upper panel) supports the predicted independence of 
A for heavier targets. They do not support the an- 
ticipated enhancement of the asymmetries compared to 
the free proton asymmetries for spin-0 and spin-1/2 nu- 
clei. However, the measured amplitude for the coherent- 
enriched sample receives contributions from incoherent 
scattering which is expected to diminish i?LU ■ The value 
^?lu = 0.91 ± 0.19 for the coherent-enriched samples 
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FIG. 6: Nuclear-mass dependence of the cos <j> amplitude of 
the beam-charge asymmetry for the coherent-enriched (upper 
panel) and incoherent -enriched (lower panel) data samples for 
hydrogen, krypton and xenon. The coherent-enriched sam- 
ples have a purity of about 67% and the incoherent-enriched 
samples a purity of about 60%. The inner error bars represent 
the statistical uncertainty and the full bars the quadratic sum 
of statistical and systematic uncertainties. 



should therefore be compared to a prediction involving 
a mixture of asymmetry amplitudes for coherent and in- 
coherent processes. For an average purity of 67% for 
the coherent-enriched samples of nitrogen to xenon (see 
Table |TTJ) and assuming that the asymmetry from the in- 
coherent portion of the yield is the same as for hydrogen, 
the predicted ratio Rlxj = 5/3 for the pure coherent pro- 
cess becomes 1.45. In Ref. [15j, both coherent and inco- 
herent scattering have been considered and the predicted 
ratio i?LU of about 1.8 for the pure coherent process be- 
comes about 1.6. These values exceed the measured ratio 
by more than three standard deviations of the total ex- 
perimental uncertainty. As shown in Figs. [5] and for 
coherent-enriched samples both beam-charge and beam- 
hclicity amplitudes for hydrogen, krypton and xenon are 
essentially independent of t within uncertainties. 

In Ref. [23[, mesonic degrees of freedom were also 
considered in the description of coherent scattering on 
nuclei and in the explanation of the generalized EMC 
effect. Such a contribution is predicted to significantly 
enhance the real part of the DVCS amplitude, which 
translates into a strong nuclear-mass dependence of 
the beam-charge asymmetry. In the absence of meson 
exchange, this asymmetry is expected to be essentially 
independent of A for heavier nuclei. The nuclear 
beam-charge amplitudes shown in Fig. [6] (upper panel) 
do not show any enhancement about the amplitude for 
the free proton and do not exhibit any dependence on A. 

In summary, the nuclear-mass dependence of az- 
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FIG. 7: Nuclear-mass dependence of the sm<f> amplitude of 
the beam-helicity asymmetry for the coherent-enriched (up- 
per panel) and incoherent-enriched (lower panel) data sam- 
ples. See Fig. [5] for the meaning of open and full circles. The 
coherent-enriched samples have a purity of about 67% ex- 
cept for He with 34%, and the incoherent-enriched samples 
a purity of about 60%. The inner error bars represent the 
statistical uncertainty and the full bars the quadratic sum 
of statistical and systematic uncertainties. This amplitude 
is subject to an additional 3.4% maximal scale uncertainty 
arising from beam polarimetry. 



imuthal beam-helicity asymmetries in electroproduction 
of real photons is measured for the first time for targets 
ranging from hydrogen to xenon. For hydrogen, krypton 
and xenon, data were taken with both beam charges 
and beam helicities allowing a separation of the shi</> 
amplitude of the squared DVCS and the interference 
terms. Also, the cost/) amplitude of the beam-charge 
asymmetry has been evaluated for those targets. This 
amplitude is consistent with earlier measurements for 
hydrogen 35J, |45j. For the coherent-enriched data 
amplitude is found to be consistent 



sample, the cos ( 
with zero for all nuclear targets, while it amounts to 0.1 
for the incoherent-enriched data sample, in either case 
not exhibiting any dependence on A within experimental 
uncertainties. 

The sin0 amplitude of the beam-helicity asymmetry 
sensitive to the squared DVCS amplitude is consistent 
with zero for all targets. The sin</> amplitude of the 
beam-helicity asymmetry sensitive to the interference 
term is significantly non-zero with a value of about —0.2 
for both the coherent and incoherent-enriched samples 
without showing any dependence on A within uncertain- 
ties. These amplitudes are compared to those of a free 
proton. The ratio i?LU = ^lu (I +) a/^lui h ^ s f° un d 
to be 0.91 ± 0.19 for the coherent-enriched sample and 
0.93 ± 0.23 for the incoherent-enriched sample, both of 
which are compatible with unity. 
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inner error bars represent the statistical uncertainty and the 
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samples is indicated for the two t bins. 
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FIG. 9: The sin0 amplitude of the beam-helicity asymme- 
try sensitive to the interference term for hydrogen, krypton 
and xenon as function of t. This amplitude is subject to an 
additional 3.4% maximal scale uncertainty arising from beam 
polarimetry. Otherwise as Fig. 8. 



No nuclear-mass dependence of the beam-charge 
and beam-helicity asymmetries is observed for heavier 
nuclei, in agreement with the general feature of models 
that approximate nuclear GPDs by a sum of nucleon 
GPDs convoluted with the distribution of nucleons in 
the nucleus. The data do not support the enhancement 
of nuclear asymmetries compared to the free proton 
asymmetries for coherent scattering on spin-0 and spin- 
1/2 nuclei as anticipated by various models 14. TBI 23j. 
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